Magnetic phase diagram of the spin-1/2 antiferromagnetic zigzag ladder 
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We study the one-dimensional spin-1/2 Heisenberg model with antiferromagnetic nearest-neighbor 
Ji and next-nearest-neighbor J2 exchange couplings in magnetic field h. With varying dimensionless 
parameters J2/J1 and h/Ji, the ground state of the model exhibits several phases including three 
gapped phases (dimer, 1/3- magnetization plateau, and fully polarized phases) and four types of 
gapless Tomonaga-Luttinger liquid (TLL) phases which we dub TLL1, TLL2, spin-density-wave 
(SDW2), and vector chiral phases. From extensive numerical calculations using the density-matrix 
renormalization-group method, we investigate various (multiple-) spin correlation functions in detail, 
and determine dominant and subleading correlations in each phase. For the one-component TLLs, 
i.e., the TLL1, SDW2, and vector chiral phases, we fit the numerically obtained correlation functions 
to those calculated from effective low-energy theories of TLLs, and find good agreement between 
them. The low-energy theory for each critical TLL phase is thus identified, together with TLL 
parameters which control the exponents of power-law decaying correlation functions. For the TLL2 
phase, we develop an effective low-energy theory of two-component TLL consisting of two free bosons 
(central charge c = 1 + 1), which explains numerical results of entanglement entropy and Friedel 
oscillations of local magnetization. Implications of our results to possible magnetic phase transitions 
in real quasi-one-dimensional compounds are also discussed. 
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PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Cx 



I. INTRODUCTION 



Frustrated quantum antifcrromagnets have long been 
a subject of active research, since Anderson^ suggested 
resonating-valence-bond ground state for a triangular 
lattice antifcrromagnct. Recent experimental studies 
of quasi-two-dimensional compounds, such as the or- 
ganic Mott insulator^ k-(BEDT-TTF) 2 Cu 2 (CN) 3 and 
the transition metal chloride Cs2CuCi4^ have further 
prompted theoretical research of anisotropic triangular 
lattice antifcrromagnets ^ri 1 : In these anisotropic quasi- 
two-dimensional antifcrromagnets combination of frus- 
trated exchange interactions and strong quantum fluc- 
tuations suppresses tendency toward conventional mag- 
netic orders, thereby opening up possibilities of exotic 
quantum states. 

A zigzag spin ladder is a one-dimensional (ID) strip of 
the anisotropic triangular lattice spin system, and can be 
regarded as a minimal, toy model of (strongly anisotropic 
quasi-two-dimensional) frustrated quantum magnets. 
Furthermore, the ID J1-J2 Heisenberg model on the 
zigzag ladder is in itself a good model for various quasi- 
1D magnetic compounds, such as (N2H5)CuCl3^ r — 
Rb2Cu 2 Mo 3 Oi2^ and LiCuVO^"-^ Despite its sim- 
plicity, the ID J1-J2 Heisenberg model has been shown 
to exhibit various unconventional phases under magnetic 
field (as we summarize below) ,2ir— In this paper we 
aim to clarify the nature of the phases in the ground- 
state phase diagram of the ID spin-1/2 J1-J2 Heisen- 
berg model under magnetic field, when both the nearest- 
and next-nearest-neighbor exchange couplings are anti- 
ferromagnetic (AF). To this end, we study in detail spin 



correlations in each phase using the numerical density 
matrix renormalization group (DMRG) method as well 
as low-energy effective theory based on bosonization. 

The Hamiltonian of the J1-J2 Heisenberg zigzag spin 
ladder is given by 
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where s/ is a spin-1/2 operator at lih site, J\ and J 2 
are respectively the nearest- and next-nearest-neighbor 
exchange couplings ( J\ > and J2 > 0), and h is external 
magnetic field along the z-direction. 

In the classical limit, the ground state of the zigzag lad- 
der J1-J2 Heisenberg antifcrromagnct has a helical mag- 
netic structure 
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and a canting angle 
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for J2/J1 > 1/4, whereas the ground state has canted 
Neel order for J2/J1 < 1/4. 

In the quantum (s = 1/2) case, the ground-state 
properties of the model ([1]) change drastically from the 
classical spin state. The ground state at zero mag- 
netic field h = has been understood quite well. For 
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small J2/J1 < {'h/Jijc-, the ground state is in a criti- 
cal Tomonaga-Luttinger liquid (TLL) phase with gapless 
excitations. The model undergoes a quantum phase tran- 
sition at (J 2 /Ji) c = 0.2411)2^— to a gapped phase with 
spontaneous dimerization 2 ^— for J 2 /J\ > {J 2 /Ji) c - It 
is also known that the model exhibits a long-range or- 
der (LRO) of vector chirality in the case of anisotropic 
exchange couplings £2r— 

With applied magnetic field, the phase diagram be- 
comes even richer. From numerical studies of the mag- 
netization process, it has been found that for a cer- 
tain range of J 2 /J\ the magnetization curve exhibits a 
plateau at one-third of the saturated magnetization and 
cusp singularities ! 21 ' 35 " — In this 1/3-plateau phase, the 
ground state has a magnetic LRO of up-up-down struc- 
ture. Furthermore, it was found that away from the 
1/3-plateau and at J 2 / J\ > 1, the total magnetization 
5 t 2 ot = J2i s i changes in units of AS* ot = 2, indicating 
that two spins form a bound pair and flip simultaneously 
as the field h increases! 21 ' 36 These characteristic changes 
in the magnetization process give accurate estimates of 
phase boundaries, which divide the parameter space into 
several regions (see Fig. [T] below), although the magne- 
tization process alone cannot give much information on 
the nature of each phase. 

Another interesting feature of the Ji-J 2 zigzag ladder 
in magnetic field is a field-induced LRO of the vector 
chirality, 

K; (n) = (a, x s l+n f. (5) 

In zero field, the vector chiral LRO has been found 
when and only when the system has an easy-plane 
anisotropy< 3 - 2 - ~ 34 i 38 i 39 In this case, due to the anisotropy, 
symmetry of the system in spin space is lowered from 
isotropic SU{2) to U(l) x Z 2 , where the U(l) and Z 2 
symmetries correspond to the rotation in the easy plane 
and the sign of pitch angle of helical spin order, respec- 
tively. While the continuous U(l) symmetry is preserved 
in the quantum case s = \s\ < 00^ the discrete Z 2 sym- 
metry can be spontaneously broken even in the quan- 
tum limit s = 1/2, thereby resulting in the vector chiral 
phase. This line of symmetry consideration suggests that 
the magnetic field, which induces the same symmetry re- 
duction, should also lead to the spontaneous symmetry 
breaking of the Z 2 symmetry. Indeed, this possibility 
was first pointed out by Kolezhuk and Vekua^ who have 
predicted from a field-theoretical analysis that the vector 
chiral LRO may set in for a large J 2 j J\ regime. Recently, 
the appearance of the vector chiral LRO under magnetic 
field was verified numericall y 42 ' 43 

In this paper, we report our numerical and analytic re- 
sults of the ground-state properties in the various phases 
that appear under magnetic field. From a thorough 
comparison of long-distance behavior of correlation func- 
tions, we identify effective theories that describe the low- 
energy physics of each phase. For this purpose, we cal- 
culate numerically various correlation functions, which 
include longitudinal-spin, transverse-spin, vector chiral, 



and nematic (two-magnon) correlation functions using 
the DMRG methodj 44 ' 45 Comparing the numerical re- 
sults with asymptotic forms derived from bosonization 
analysis, we find that, in addition to the gapped dimcr 
phase, 1/3-plateau phase, and fully polarized phase, the 
system exhibits four critical phases: (i) a phase with 
one-component TLL which is adiabatically connected to 
the ground state of the ID Hcisenbcrg antifcrromagnet 
(TLL1 phase), (ii) a two-component TLL phase (TLL2 
phase), (iii) a vector chiral phase, and (iv) a spin-density- 
wave phase with two-spin bound pairs (SDW2 phase). 
The low-energy states in the TLL1, vector chiral, and 
SDW2 phases turn out to be one-component TLLs (a con- 
formal field theory with central charge c = 1). Further- 
more, we provide quantitative estimates of non-universal 
parameters appearing in the low-energy effective theo- 
ries, such as the TLL parameter and incommensurate 
wavenumbers of spin correlations, as functions of J 2 / J\ 
and the magnetization. In particular, our results of the 
TLL parameter, which controls decay exponents of cor- 
relation functions, have direct relevance to experimental 
observables, e.g., a magnetic LRO emerging in real quasi- 
1D compounds with weak interladdcr couplings and tem- 
perature dependence of relaxation rates (1/Ti) in nuclear 
magnetic resonance experiments 1 46 1 47 We also propose a 
two-component TLL theory to describe the TLL2 phase. 

This paper is organized as follows: In Sec. UH we show 
the ground-state phase diagram under magnetic field (sec 
Fig. P, which contains the TLL1, 1/3-plateau, SDW 2 , 
vector chiral, TLL2, dimer, and fully polarized phases. 
We briefly summarize the characteristics of each phase. 
In the following sections, we discuss in detail our numeri- 
cal results for correlation functions and effective theories 
for each phase. In Sec. IIII1 we consider the TLL1 phase, 
which appears in small J 2 /J\ regime. The correlation 
functions obtained with the DMRG method are shown to 
be fitted well to analytic forms obtained from a bosoniza- 
tion theory for a weakly-perturbed single Heisenberg spin 
chain, and the decay exponents of the spin correlation 
functions are estimated accurately. This analysis reveals 
that the dominant correlation function changes from the 
staggered transverse-spin correlation to incommensurate 
longitudinal-spin one as J 2 /J\ increases. In Sec. IIV1 we 
discuss the SDW2 phase, which appears at larger J 2 j J\. 
From the fitting of numerical data to bosonization theory, 
we show that the low-energy excitations are described 
by a one-component TLL with quasi-long-ranged domi- 
nant incommensurate longitudinal-spin and subleading 
nematic correlations and short-ranged transverse-spin 
correlation. Section [V] discusses the 1/3-platcau phase. 
We show that the numerically found up-up-down spin 
structure is understood in terms of the bosonization the- 
ories for the neighboring TLL1 and SDW2 phases. In Sec. 
IVI1 we consider the vector chiral phase, which is also a 
one-component TLL. The fitting analysis shows that the 
vector chiral phase is characterized by the vector chiral 
LRO and the incommensurate quasi-LRO of the trans- 
verse spins. In Sec. I VIII we develop a two-component 
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TLL theory, i.e., two free boson theories (central charge 
c = 1 + 1), as a low-energy effective theory for the TLL2 
phase. We confirm the central charge c = 2 through 
numerical computation of entanglement entropy. The 
consistency between the effective theory and the DMRG 
result is shown by examining a few dominant Fourier 
components in the local magnetization profile near open 
boundaries. Section IVIIII contains summary and discus- 
sions on implications of our results to real quasi-lD com- 
pounds with weak interladder couplings. 

II. PHASE DIAGRAM 

Figure[T]presents the magnetic phase diagram based on 
the numerical results obtained in this paper as well as in 
previous studies. The diagram is shown in the Jif J\ ver- 
sus magnetization M plane in Fig. HJa) and in the J<ij J\ 
versus h/J\ plane in Fig. QJb), where M = (l/£)X)z s i 
is the magnetization per site and L the system size. The 
system exhibits at least four critical phases, i.e., TLL1, 
TLL2, vector chiral, and SDW2 phases, in addition to 
three gapped phases including the dimer phase at M = 0, 
the 1/3-plateau phase (M = 1/6), and the fully polarized 
phase (M =1/2). 

It has been revealed that the magnetization process 
of the zigzag ladder ([T]) has remarkable features; 21 ' 35 " — 
for small J2/J1 (0.25 < J2/J1 < 0.7), the magnetiza- 
tion curve has at most two cusp singularities at higher 
and lower fields, h = h c \ and h c2 , which correspond 
to boundaries between the TLL1 and TLL2 phases. A 
magnetization plateau also appears at M = 1/6 for 
0.487 < J2/J1 < 1.25 and h pl < h < h p2 2LZL For 
large J2/ J\ , the magnetization process exhibits two-spin 
flips with AS^ ot = 2 in an intermediate field region 
h m i < h < /i m2 . See Figs. 2 and 3 in Ref.[2l]for these re- 
sults. At zero magnetization, the ground state is gapless 
for J2/J1 < 0.2411 and dimerized for 0.2411 < J 2 /Ji- 
The spin gap in the dimerized phase vanishes at a crit- 
ical field /id- The ground state is fully polarized above 
the saturation field h s . The critical fields h c ±, h c2 , h p \, 
hp2i h m \, /i m 2 7 ^d, and h s are plotted in Fig. [TJb) with 
solid lines. 

To reveal the nature of ground states in each region, 
we have calculated several correlation functions, using 
the DMRG method, for the system with up to L = 160 
spins with open boundaries. We have kept typically 300 
block states in the calculation (up to 400 states for some 
cases), and confirmed the convergence of the calculation 
by checking the dependence of results on the number of 
kept states. We have calculated the longitudinal-spin 
correlation function (s^sf,), the transverse-spin correla- 
tion function (sfsf,), the vector chiral correlation func- 
tion (aw «;/ ) with n,n' — 1,2, the nematic (two- 
magnon) correlation function (s^~ s^ +1 s^ sf, +1 ) , and the 
local spin polarization (sf), where (•■■) denotes the ex- 
pectation value in the ground state. To lessen the open- 
boundary effects, we have computed the two-point corre- 
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FIG. 1: (Color online) Magnetic phase diagram of the spin- 
1/2 antiferromagnetic zigzag ladder (a) in the J2/J1 versus 
M plane and (b) in the J2/J1 versus h/Ji plane. In (a), sym- 
bols represent parameter points for which their ground-state 
phases are identified: Open (o) and solid (•) circles represent 
the TLL1 phase with dominant transverse- and longitudinal- 
spin correlation, respectively. Diamonds (o), triangles (A), 
and squares (□) respectively represent the SDW2, vector chi- 
ral, and TLL2 phases. The solid line shows the 1/3-plateau 
phase. The dotted curves are the guide for the eye. In (b), 
symbols P, VC, D, and F indicate the 1/3 plateau, vector chi- 
ral, dimer, and fully polarized phases, respectively. T he p hase 
boundaries shown by solid lines are obtained in Ref. [2l| from 
the numerical results of magnetization curves, except for the 
boundaries between the vector chiral and TLL2 phases which 
are obtained from the analysis of correlation functions in the 
present paper. The dotted line in the TLL1 phase represents 
the crossover line between the transverse- and longitudinal- 
spin dominant regimes. 



lation functions for several pairs of (1,1') with fixed dis- 
tance r — \l — l'\ and taken their average for the estimate 
of the correlation at the distance r. In the following, 
we use the notation (■ • • ) av for the averaged correlation 
functions. 

Figure [2] shows typical spatial dependence of averaged 
correlation functions in the critical phases. We note 
that the bending-down behaviors of the averaged cor- 
relation functions seen for large distance (e.g., r > 100 
for L = 160) are due to boundary effects and should not 
be confused with intrinsic behaviors in the bulk. Ana- 
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FIG. 2: (Color online) Typical spatial dependence of correlation functions in critical phases; (a) TLL1 phase (transverse-spin 
correlation dominant), (b) TLL1 phase (longitudinal-spin correlation dominant), (c) SDW2 phase, (d) vector chiral phase, and 
(e) TLL2 phase. Absolute values of the averaged correlation functions are plotted. 



lyzing the long-distance behavior of correlation functions 
in each parameter regime, we have determined the low- 
energy effective theory for each phase. The parameter 
points in the phase diagram at which numerical results 
are explained successfully by the effective low-energy the- 
ory of the corresponding phase are shown with symbols in 
Fig.fTJa). We summarize properties of each phase below. 

TLL1 phase: In small J2/ J\ regime, the ground state is 
adiabatically connected to the one-component TLL of the 
antifcrromagnetic Heisenbcrg chain with only J\ under 
magnetic field. For relatively large J 2 / ' J\ (0.25 < J2/J1) 
the boundaries of the TLL1 phase are defined by the 
cusp singularities in the magnetization curvei 21 i 35 In this 
phase, both the longitudinal-spin fluctuation (sqS*) — 
(sq}(s*) and transverse-spin correlation functions (s'qS*) 
decay algebraically. The former shows incommensurate 
oscillations with a wavenumber Q = n(l ± 2M), while 
the latter is staggered, Q = n. The numerical estima- 
tion of the decay exponents, shown in Sec. IIII1 indicates 
that the dominant correlation function changes from the 
staggered transverse-spin correlation to incommensurate 
longitudinal-spin correlation as J2/J1 increases [sec Figs. 
[2] (a) and (b)]. The TLLI phase is thus divided by the 
crossover line into two regions of different dominant cor- 
relations, as shown in Fig. [TJ 

SDW2 phase: For large J2/J1, there is a phase 
where the magnetization process changes by the steps 
of AS* t z ot = 2SL We show in Sec. |IV] that this phase is 
described by a one-component TLL theory, which was 
originally derived from the weakly-coupled AF Heisen- 
berg chains in the limit J 2 /Ji » h 22 ' 41 ^^ The phase 
is characterized by the quasi-long-ranged longitudinal- 



spin and ncmatic correlation functions, (sqS^ 



and (sq Si s~ s~ +1 ) , which are dual to each other, and 
by the short-ranged transverse-spin correlation function 
(sgS^) reflecting a finite energy gap to single-spin- flip 
excitations, as shown in Fig. [5{c). The longitudinal- 
spin correlation is incommensurate with the wavenum- 
ber Q2 = ±7r(l/2 + M). Numerical analyses of corre- 
lation functions reveal that the longitudinal-spin corre- 
lation function is dominant in the whole parameter re- 
gion of this phase. We thus call this phase the SDW2 
phase. We note that the same phase has been found in 
the zigzag ladder (JXJ) with ferromagnetic J\ and AF J 2 
as well, 2 ^^ 3 . 

1/3-plateau phase: At one third of the saturated mag- 
netization, M — 1/6, there is a magnetization-plateau 
phase in the intermediate parameter region 0.487 < 
J2/J1 < l-25 i 21 ' 37 This phase is characterized by a field- 
induced excitation gap and a spontaneous breaking of 
translational symmetry accompanied by a magnetic LRO 
of the up-up-down structure. The ground state is three- 
fold degenerate. As shown in Sec. [V] all two-point corre- 
lation functions exhibit exponential decay, in accordance 
with the fully-gapped nature of the phase. 

Vector Chiral phase: The vector chiral phase is char- 
acterized by the LRO of the vector chirality as well 
as quasi-LRO of incommensurate transverse spins, which 
decays algebraically in space. The discrete Z2 symme- 
try corresponding to the parity about a bond center is 
broken spontaneously and the ground state is doubly de- 
generate in the thermodynamic limit. This vector chi- 
ral state is a quantum counterpart of the classical heli- 
cal state. Though the classical helical state appears in 
1/4 < J2/J1 for arbitrary magnetization, the quantum 
vector chiral phase is found only in two narrow regions 
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separated by the SDW2 and 1/3-platcau phases ; ' see 
Fig. [T] We show that the vector chiral phase is also 
described by a one-component TLL theory which can 
be formulated starting from the two weakly-coupled AF 
Heisenberg chains for J2/J1 S> li 32 ' 41 The correlation 
functions in this phase will be discussed in Sec. IVII 

TLL2 phase: The TLL2 phase occupies two parame- 
ter regions adjacent to the TLL1 phase and the vector 
chiral phase. The TLL2 phase is described as two Gaus- 
sian conformal field theories (central charge c = 1 + 1), 
or a two-component TLL, having two flavors of free 
massless bosonic fields as its low-energy excitations. In 
the Jordan- Wigner fermion representation, fermions have 
two separate Fermi seas, and the two bosonic fields rep- 
resent particle-hole excitations near the two sets of Fermi 
points. In the TLL2 phase all correlation functions de- 
cay algebraically and have incommensurate wave num- 
bers which are linear functions of the two Fermi momenta 
of Jordan- Wigner fermions. We will discuss these prop- 
erties and the low-energy effective theory in Sec. IVIII 

Dimer phase: For J 2 /Ji > 0.2411 and at M = 0, 
the ground state of the J1-J2 AF Heisenberg zigzag spin 
ladder is spontaneously dimerized<^r— The ground state 
is doubly degenerate in the thermodynamic limit, and 
there is a gap to lowest excitation. 

Fully polarized phase: When applied magnetic field is 
larger than the saturation field, h > h s , the ground state 
is in the fully polarized phase with saturated magnetiza- 
tion M = 1/2. As the field decreases, the fully-polarized 
ground state is destabilized by softening of single- magnon 
excitations, which have the dispersion, 



£k = Jl(cOsfc — 1) + J 2 (cOs2k — 1 



) + h. 



(6) 



When J2/J1 < 1/4, the magnon dispersion has a single 
minimum at k = Tt, while, when Jij J\ > 1/4, there 
are two energy minima at k = ± arccos(— Ji/4 J2). The 
saturation field h s is given by h s /Ji =2 for J2/J1 < 1/4 
and h s /Jx = 2J 2 /Ji + 1 + Ji/(8J 2 ) for J 2 /Ji > 1/4. 



III. TLL1 PHASE 

In this section, we discuss the TLL1 phase appearing 
for small J%j J\. Since the parameter space of this phase 
includes the AF Heisenberg chain with 3 2 = 0, we natu- 
rally expect that the TLL1 phase should share the same 
properties with the single Heisenberg chain. Here, we 
first briefly review the TLL theory for the AF Heisen- 
berg zigzag ladder with weak J2 coupling. We then com- 
pare the theory with the numerical results of correlation 
functions for the zigzag ladder (fTJ) with J 2 > 0. 

It is well known that the low-energy properties of a sin- 
gle Heisenberg chain under magnetic field (|M| < i and 
J2 = 0) is described as a TLL.— Since the (leading) op- 
erator generated from weak J 2 coupling is irrelevant in 
applied magnetic field (and marginally irrelevant with- 
out magnetic field) in the renormalization-group sense, 



the low-energy effective theory for small J2/J1 is adia- 
batically connected to the TLL theory of the single AF 
Heisenberg chain (J2 = 0). Hence the low-energy exci- 
tations in the TLL1 phase are free massless bosons gov- 
erned by the Gaussian model, 



Ho 



2 ' 



dx) K V dx 



(7) 



where (</>, 9) are bosonic fields satisfying the equal-time 
commutation relation [<p{x),d y 6(y)\ = iS(x — y). The 
TLL parameter if is a function of J2/J1 and M. We 
have taken the lattice spacing to be one and identify the 
continuous coordinate x with the site index I. The spin 
velocity v is of order J\, except for the saturation limit 
M — > 1/2, where v —¥ 0. The spin operators s/ can be 
expressed in terms of the bosonic fields as 

1 d<f>(x) 



Sf=M; 

V7r dx 

- (-l)Wn[27rMZ + V4n<f>(x)] + • • • , (8) 
s+ = {-l) l be 1 ^ 9 ^ 

+ b'e l ^ e{x) sin[27rA/Z + V^<f>(x)] + ■■■ , (9) 

where a, b, and b' are nonuniversal positive constants, 
whose numerical values are known at J 2 = 0.— r — Equa- 
tions (|7|), ©, and © define the effective theory for the 
TLL1 phase, with which asymptotic forms of spin corre- 
lation functions arc obtained as 



(t&st) = M 2 — 71 



47T 2 7 



_(-l) r cos(27rMr) , . 



\r\v 



s o r 



:(-l) r i4? _cos(27rMr) 



|r|!A; 



\r\ri+Vv 



(11) 



where A\ = a 2 /2, A% = b 2 /2, Af = b' 2 j^ (with appro- 
priate short-distance regularization), and the decay ex- 
ponent r\ is related to the TLL parameter K by rj = 2K. 
Equations (jTU)) and (fTTj) tell us that for 77 > 1 the stag- 
gered transverse-spin correlation function (sfsf,) is dom- 
inant, while the incommensurate longitudinal-spin cor- 
relation (sfsf,) with a wavenumber Q = tt(1 ± 2M) is 
dominant for 77 < 1. At J 2 = the decay exponent r\ can 
be calculated exactly using Bethe ansatzj 57 i 58 r\ increases 
monotonically as M increases, from 77 = 1 at M = to 
i] = 2 for M — > 1/2. Therefore, at J 2 = 0, the transverse- 
spin correlation (sqS*) is always the most-slowly decay- 
ing one for < M < 1/2. For finite J2 > 0, the exact 
value of the exponent 77 is known in the limit M — ¥ 0. For 
J2/J1 < 0.2411 where the ground state at M = is in the 
TLL1 phase, r\ = 1 at M = because of the SU(2) sym- 
metry. On the other hand, for J2/J1 > 0.2411, i.e., when 
the ground state at M = is in the dimer phase ; 59 ' 60 
i] —> 1/2 as M —> 0. This means that 77 is singular at 
J2/J1 = (J 2 /Ji) c and M = 0. 

One can also derive the (same) effective theory for 
the TLL1 phase, starting from the saturation limit for 
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2ttLM 
L + l-25' 



and 



21L + 1-26) 
sin 



2(L + 1 - 26) 



< J2/J1 < 1/4. In this limit the system can be viewed where 
as a dilute gas of interacting hard-core bosons (magnons) 
with one flavor, as the magnon dispersion ((6]) has a sin- 
gle minimum at k = n. The hydrodynamic theory for 
the one-flavor interacting bosons is nothing but the TLL 
theory, Eq. ([7]).— This approach naturally gives the same 
asymptotic forms of spin correlators as Eqs. © and 
Furthermore, in the saturation limit M — > 1/2, 77 — > 2 in 
the TLL1 phase (i.e., J2/J1 < 1/4), since the dilute limit 
of the hard-core bose gas is equivalent to a free fermion 
gas. f a (x) 

Next we discuss our DMRG results of the transverse 
and longitudinal spin correlation functions (sfsf,) and 

(sfsf,) and the local spin polarization (sf). To achieve Jn ^ ^ , <<; L ^ Eq ™j reduces tQ 
better numerical convergence and efficiency, the DMRG 
calculation was done for finite systems (L spins) with 
open boundaries. We thus compare the numerical re- 
sults with the correlation functions calculated analyti- 
cally from the effective theory ([7]) by imposing appropri- 
ate boundary conditions on the bosonic field (f>^r— To 
this end, we have taken the Dirichlct boundary conditions 
(f)(6) = (f>(L + 1 — 6) = Of£i where <5 is a free parameter to 
be determined later. For example, spatial dependence of 
the magnetization is given by 



(13) 



(14) 



M 



(-l)'o 
[2(1 - 6)]n/ 2 



sin[2vrM(Z - 6)} 



(15) 



W/ ^ ' 2tt f v/2 (2(l-6)) 



(12) 



The presence of open boundaries gives rise to "Friedel 
oscillations" in the local magnetization. The wave num- 
ber of the oscillations is "2fcp" of the Jordan- Wigncr 
fcrmions, which equals Q = ir(l ± 2M) for L ^> 1. Sim- 
ilarly, the longitudinal and transverse spin correlation 
functions are modified by boundary contributions as 



= Z(l,l';q) 



4tt 2 



h(i-v) 
(-if 



h(l + V-26) 



qa 
2^ 



(-1) 1 sm[q(l - 6)} (-1) 1 ' wo[q(l> - 6)} 



f v/2 {W-s)) 



M-V a 2 



2f v/2 (2(l-6))f v/2 (2(l'-6)) 

rjaj (-1) 1 cos[q(l - 6)} 
2tt 



cos [ q (l _ l')]lA±L-^L _ C0B [ q (l + 26)} 



f v/2 (2(U - 6)) 

Ml -i') 



Mi + 1' -26) 



[g(l + 1'- 26) + g(l - I')} + 



(-1) 1 ' cos[q(l'-6)] 



/„/a(2(J-*)) 
(st)(st,) = Z(l,l';q)-z(l;q)z(l , ;q), 
X(l,l';q) 

f 1/2 r 1 (2(l-6))f 1/2 , 1 (2(V-6)) j (-l)'- l 'b 2 S gn(l-l')bb? 

h/ v (i-i')h /v (i + i'-^) \ 2 2 
b' 2 



f n/2 (2(U - 6)) 



[g(l + 1'- 26) - g(l - l')\ 



(16) 
(17) 



(-1)' cos{q(l-8)\ _ {-Vfco^^-S)\ 
/„/ 2 (2(Z-<S)) fn/ziW-S)) 



Af v/2 (2(l-6))f v/2 (2(l'-6)) 



cos[q(l + I' - 26)} 



Mi -i') 

Ml + V - 26) 



cos[q(l — I') 



Mi + 1' -26) 



mi - n 



(18) 



r 



where 



2(L + 1 - 26) 



cot 



nx 



2(L + l-26) 



(19) 



In the limit \L/2 - l\ < L and \L/2 - l'\ < L, boundary 
effects go away, and Eqs. (jTHJ) and (JT5J reduce to Eqs. 
(JTUJ) and (fTTj) . In the fitting procedure discussed below, 



we have optimized 6 to achieve the best fitting of (sf) 
and (sfsf) — (sf)(sf), whereas we set 6 = for (sfsf,) as 
it has turned out that the numerical data of (sfsf,) can 
be fitted sufficiently well without optimizing 6. 

Figured shows DMRG data of (sf), (sfsf) - (sf)(sf), 
and (sfsf) for (J 2 /J^M) = (0.1,0.375) and (0.5,0.125). 
In the same figures, we show the fits to Eqs. (fT5|) . (|17l) . 
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FIG. 3: (Color online) Correlation functions in the antiferromagnetic zigzag ladder with L — 160 spins in the TLL1 phase; 
(a) local spin polarization (sf), (b) longitudinal-spin fluctuation (sfsf/) — (sf)(sf/), and (c) transverse-spin correlation function 
(sfsf). The upper and lower panels show the results for (J2/ Ji, M) = (0.1,0.375) and (0.5,0.125), respectively. The open 
symbols represent the DMRG data and the solid lines and circles are fits to Eqs. <|12[) . (|17|) . and (JTHJ . In (b) and (c), the data 
for l = L/2- [r/2] and I' = L/2 + [(r + l)/2] are shown as a function of r = |Z - Z'|. 



and (|18jl . Clearly, the fits are in excellent agreement with 
the numerical results. We emphasize that only three fit- 
ting parameters, 77, 5, and a (77, b, and 6') are used in the 
fitting of (sf) and (sfsf,)— (sf)(sf,) {(sfsf,)). Wehaveob- 
taincd almost the same good quality of fits for the param- 
eter points marked by open and solid circles in Fig.[TJb), 
which cover almost the entire region of the TLL1 phase. 
These results thus demonstrate that the TLL1 phase is 
described by the effective TLL theory given by Eqs. (J7J, 
(|HJ), and which is indeed the same TLL theory as 
that of the AF Heisenberg chain (J 2 = 0). 

Figure [4] shows dependence of the exponent 7/ on the 
magnetization M in the TLL1 phase, obtained from the 
fitting of the transverse-spin correlation (sfsf,). Similar 
estimates of 77 are obtained from the other correlators 
(not shown) . For small J2/J1, 77 exhibits essentially the 
same behavior as a function of M as 77 (Af) at J 2 = 0; for 
J2/J1 i5 0.15, 77 increases monotonically from the univer- 
sal value 77 = 1 at M = to 77 = 2 at M -> 1/2 as M 
increases. In this regime the transverse-spin correlation 
(sgsf) is dominant for any M. The situation changes 
as J2/J1 gets larger. With increasing J2/J1, 77 decreases 
and becomes smaller than 1 at J2/J1 = 0.2 for interme- 
diate magnetization M. As J2/J1 is further increased 
in the TLL1 phase, the exponent 77 gets smaller than 
1 for any M. Thus, the system undergoes a crossover 
from the small J2 / Ji region with the dominant stag- 
gered transverse-spin correlation to the large J2 / Ji re- 
gion where the incommensurate longitudinal-spin corre- 
lation with Q = 7r(l ± 2M) is dominant. The crossover 
line is shown in the phase diagram, Fig. Q] The result is 
consistent with the earlier study^ in which 77 was esti- 
mated at M = 1/6, 1/4 and 1/3 for small systems. Such 
a crossover between ground states with the different dom- 



inant spin correlations has also been found for the J1-J2 
zigzag ladder with bond alternation^ - — 

As mentioned above, 77 is expected to approach 1/2 as 
M -> for J2/J1 > (J 2 /Ji) c = 0.2411. Our numerical 
results at J2/J1 = 0.5 are consistent with this theoretical 
prediction. However, as J2/J1 approaches (h/Jijc from 
above, the value of 77 at the smallest M = 0.025 becomes 
larger toward 77 = 1, the value expected for J2/J1 < 
(J2I J\)c- This implies that 77 increases very rapidly from 
1/2 at small M for this parameter regime of J2/J1 < 0.3, 
where the spin gap in the dimer ground state at AI = 
is exponentially small (thereby small M is sufficient to 
wipe out dimer instability). 

The data points for 0.3 < J2/J1 < 0.5 end at the 
boundary to the TLL2 phase for larger M. Our results 
seem to indicate that 77 changes continuously along the 
TLL1-TLL2 phase boundary. 

When the magnetization is close to M = 1/6, the 
TLL1 phase has an instability to the 1/3-platcau phase. 
In the Jordan- Wigner fermion picture, the instability is 
caused by umklapp scattering of three fcrmions, and the 
1/3-plateau phase corresponds to a density wave state 
of the fermions i 67 ' 68 The three-particle umklapp scatter- 
ing is irrelevant at small J2/J1 but becomes relevant for 
larger J2/J1. This explains why the 1/3-plateau phase 
emerges at J2/J1 > 0.5 in the phase diagram (Fig. [!}, as 
we discuss below. 

The effective Hamiltonian yielding the 1 /3-plateau has 
the formS-es 



H = Ho + X J dxsin|7r(6M - l)x + 3V47T(/)(a;)J , (20) 

where Ho is the Gaussian model (J7J for the TLL1 phase 
and A is the coupling constant for the three-particle umk- 
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FIG. 4: (Color online) M dependence of the exponent 77 
for the TLL1 phase estimated from the fitting of (sfsfi) for 
the antiferromagnetic zigzag ladder with L = 160 spins. The 
error bars represent the difference of the estimates obtained 
from the fitting of the data of different ranges. The vertical 
dashed line corresponds to M — 1/6 where the 1/3-plateau 
can appear for large J2/J1. The dotted line at 77 — 1 rep- 
resents the boundary between the regimes of dominant stag- 
gered transverse-spin correlation (77 > 1) and dominant in- 
commensurate longitudinal-spin correlation (77 < 1); see Eqs. 
(|10[) and (|11[) . The exponent 77 relates to the parameter K as 
77 = 2K in the TLL theory for the TLL1 phase. 



lapp scattering. The umklapp term is accompanied by 
an oscillating factor with a wavenumber it(6M — 1) = 
37r(l + 2M) and becomes uniform at M = 1/6. If 
we fix the magnetization at M = 1/6 and increase 
J2/J1, then the three-particle umklapp term becomes 
relevant for K < 2/9 (77 < 4/9). Indeed, we see in 
Fig. |4]that the estimates of 77 near M = 1/6 are larger 
than 4/9 for J2/J1 < 0.4 and become close to 4/9 at 
J2/J1 = 0.5. This result is consistent with the esti- 
mated critical value (J 2 /Ji) p i = 0.487 which was ob- 
tained from the analysis of the level spectroscopy in 
Ref . l37l. For J2/J1 > {Jz/Jiipi we can approach the 1/3- 
plateau phase by changing the magnetic field h. This is in 
the universality class of commensurate-incommensurate 
transition ! 59 ! 70 In this case we expect that, as M — >■ 1/6, 
the TLL parameter K approaches 1/9, or, cquivalently, 
T] —> 2/9 £2- On the other hand, our numerical data for 
J2/J1 = 0.6 seem to be much larger than the theoreti- 
cal value 2/9 at M — > 1/6. Although this disagreement 
might suggest that there exist rather large errors in the 
estimates of 77 for large J2/J1, we rather expect that 
77 for Jil J\ = 0.6 should actually show rapid decrease 
very close to M = 1/6 to recover the predicted behav- 
ior, 77 2/9 as M — s- 1/6. Numerical verification of this 
would require calculations on much larger systems. 



IV. SDW 2 PHASE 

In this section we discuss the SDW2 phase. This phase 
is characterized by two-spin flips AS Z = 2 in the magne- 
tization process! 21 ' 71 The parameter space of the SDW2 
phase extends to large J2/J1, see Fig. [1] Its low-energy 
effective field theory is obtained in the limit J2/J1 —> 1, 
and we will give a short review on it below i 22 i 41 ' 50 i 52 
Then, by comparing our DMRG data of correlation func- 
tions for 3\j 3-2. 1 with the analytic results, we demon- 
strate that the effective theory is valid in the whole pa- 
rameter space of the SDW2 phase, as expected from the 
principle of adiabatic continuity*^ 2 . 

In the limit J2 3> Ji , the zigzag spin ladder ([T]) can be 
viewed as two Heisenberg chains with nearest-neighbor 
exchange J2 coupled by weak interchain exchange J\ . It 
is natural to bosonize each chain separately first and then 
incorporate the interchain coupling J± perturbativcly. In 
this scheme, the original spin operators are written as 



'2j+n 



■ M 



1 d4> n (x n ) 



Fk dx 

- (-l)%sm[27rMj + V4^0„(x„)] + • • ■ , (21) 
4 + „=(-l)W^M-) 

+ b'e 1 ^ 6 ^ sin[27rMj + V^4> n {x n )] +••-, 

(22) 

where {4> n ,0n) are the bosonic fields for each chain n = 
1,2. The coordinate x is related to the site index I = 
2j + n (n = 1, 2) as x x = j - 1/4 and x 2 = j + 1/4. The 
low-energy theory of each AF Heisenberg chain has the 
same form as Ho in Eq. ([7]). The bosonized form of the 
interchain coupling J\ can be found from Eqs. (|21j) and 
(1221. We then obtain the effective Hamiltonian2 2 .>ii> 5 £>£ 2 . 



dx 



K v 



dx 



v=± 

+ gi J dx sin(v / 87T'/>- + nM) 

/d0 
dx-j±sm(V2^6-), 



1 

K v \ dx 



(23) 



where the interchain coupling gives the nonlinear inter- 
action terms with the coupling constants 

gi = .ha 2 sin(TrM), g 2 = y V^b 2 (24) 

in lowest order in J\ . Here we have introduced symmetric 
(+) and antisymmetric (— ) linear combinations of the 
bosonic fields, (f>± = (4>i ± (fa)/y/2, 6± = {6 1 ± 6» 2 )/\/2. 
In lowest order in J\ the TLL parameters K± are given 
by 



K± = K[1t 



■KV 



(25) 
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where K is the TLL parameter of the decoupled Hcisen- 
berg chains 4i This suggests that K + is less than 1 and 
decreases with J± at the limit J\ <C J 2 . The spin ve- 
locities v± are of order J 2 in the weak-coupling regime, 
except for M — > 1/2 where v± — > 0. 

The effective Hamiltonian ([23| has two competing in- 
teractions (oc g\ and g^). The fate of the ground state is 
determined by which one of the two interactions grows 
faster in renormalization-group transformations. If the 
gi term is dominant, the SDW 2 phase is realized. We 
discuss this case below. On the other hand, if the g 2 
term is most relevant, then the ground state is in the 
vector chiral phase; this case is discussed in Sec. IVIl 

Let us assume that the 171 term wins the competition. 
Then the field 4>— is pinned at a minimum of the potential 
<7i sin(v / 87r0- + 7rM), 



L=160, M=0.125 



-} 



M 



(26) 



as gi oc J 1 > 0. Since the pinned field 0_ can be taken 
as a constant, the difference of (the uniform part of) two 
neighboring spins vanishes, — s 2j+2 = \J'2/ 1T dx4>- ~ 
0. This means that the two spins are bound and explains 
the steps A5^ ot = 2 in the magnetization process^ - The 
dual field fluctuates strongly and we can therefore 
safely ignore the 172 coupling. The antisymmetric sector 
(</>_,#_) has an energy gap, which corresponds to the 
binding energy of the two-spin bound state. 

Since the bosonic fields ((/>+, #+) in the symmetric 
sector are not directly affected by the relevant inter- 
chain couplings, they remain gapless and constitute the 
one-component TLL. The effective Hamiltonian for the 
SDW 2 phase is the Gaussian model, 



v± 
2 



dx 



K 4 



d6+ 

dx 



1 



dx 



(27) 



Equations jUJ), (|35J, (gSJ), and flU} represent the TLL 
theory for the SDW2 phase. 

Straightforward calculations yield the longitudinal- 
spin and nematic (two-magnon) correlation functions in 
the thermodynamic limit, 



S 0~ S l~ S r S r+l) 



77 A\ 



cos 



\ -+M 



(28) 



(-l)^o 
| r |i/u 



cos 



?rr [ - + M 



(29) 



where the exponent rj = K + , and we have introduced 
positive numerical constants A i and A\. These corre- 
lations are quasi-long-rangcd and dual to each other. If 
X] < 1, the incommensurate SDW correlation [the third 
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FIG. 5: (Color online) Absolute values of the averaged 
transverse-spin correlation function, |(sgs^) a v|, in the antifer- 
romagnetic zigzag ladder with L — 160 spins for M — 0.125 
and several J2/J1. 



term in Eq. (|28p ] is the most dominant, while the stag- 
gered nematic correlation is the strongest for rj > 1. 
The perturbative result (f25|) indicates that the incom- 
mensurate SDW correlation is dominant (A" + < 1) for 
small Ji/J-2- We will see below that this holds true for 
J1/J2 ~ 1 as well. The wavenumber of the SDW quasi- 
LRO is Q2 = ±7r(l/2 + M), which is distinct from that of 
incommensurate correlations in other phases and is char- 
acteristic of the SDW2 phase. We note that in the SDW2 
phase of the ferromagnetic (J% < 0) J1-J2 zigzag ladder, 
the characteristic wavenumber is Q = ±tt(1/2 — M) i 22 ' 47 
Such a spin-density-wave state with the incommensurate 
wave vector is also found in the spatially-anisotropic tri- 
angular antiferromagnet in magnetic fields 

The transverse-spin correlation function (sqS^) decays 
exponentially as the operator sf includes the strongly 
disordered field. The exponential behavior is a direct 
consequence of the finite-energy cost for creating a single- 
magnon excitation and is a hallmark of the SDW2 phase. 

Let us discuss numerical results. Figure HJc) shows 
typical behaviors of the averaged correlation functions in 
the SDW2 phase. The longitudinal-spin and two-magnon 
correlation functions decay algebraically and the former 
is clearly dominant. By contrast, as shown in Fig.[5j the 
transverse-spin correlation decays exponentially. This 
can be seen as evidence for the appearance of two- 
magnon bound states in this parameter regime. The cor- 
relation length of transverse spins becomes larger with in- 
creasing J2I J\- This is in accordance with the bosoniza- 
tion prediction that the energy gap for the single-spin 
excitation is generated by the cosine term with the coef- 
ficient gi ~ J 1 for J1/J2 <C 1. We have found essentially 
the same behavior of the correlation functions as shown 
in Fig. [2jc) for the entire parameter region where the 
two-spin-flips with ASf ot = 2 are observed in the magne- 
tization process. After the dominant correlation function 
and the formation of two-magnon bound pairs, we call 
this phase the SDW 2 phase. 

In order to estimate the exponent 77 and to further 



10 




0.01 



100 / 150 
Jg/J^l.5, M=0.'l ] 




FIG. 6: (Color online) Correlation functions in the antifer- 
romagnetic zigzag ladder with L — 160 spins in the SDW2 
phase; (a) local spin polarization (sf) and (b) longitudinal- 
spin fluctuation (sfsfi) — (sf)(sf,). The upper and lower pan- 
els show the results for (J 2 /Ji,M) = (1.5,0.1) and (1.0,0.1), 
respectively. The open symbols represent the DMRG data 
and the solid symbols are the results of fitting to Eqs. (|30[) 
and (J3T). In (b), the data for I = L/2 - [r/2] and I' = 
L/2 + [(r + l)/2] are shown as a function of r = \l — 



demonstrate the validity of the effective theory for the 
SDW2 phase, we fit the DMRG data of the local-spin 
polarization (sf) and the longitudinal-spin fluctuation 



(sf)(sf,) to analytic forms obtained from the 



bosonization approach. Using Eqs. (|2T|) . (|26|) . and (|27|) 
and applying the Dirichlet boundary condition in the 
same manner as in Sec. IIII1 we obtain the correlators 



for a finite open zigzag ladder as 
(st)= 1 --2z{l;q), 



(30) 



(sfsf,) - (sf ) (sf,) = 4 [Z(l, l'; q) - z(l- q)z(l'- q)] , (31) 
with 

2ttL (\ M s 



* L + l-26 \4 2 
In the limit l<L, Eq. (J3UJ) reduces to 



= M 



2(-l)'a 
[2(1 - Sftfi* 



2 



(32) 



(1-5) , (33) 



showing Friedel oscillations with wavenumber tt(^ + M). 

Figure [5] shows DMRG results and their fits to Eqs. 
([30]) and (JSJ). The results for J 2 /J\ = 1.5 show that 
the numerical data at relatively large J2/J1 are fitted 
pretty well by the analytic forms. Note that only three 
fitting parameters, 77, a, and 8, are used in the fitting 
procedure. For smaller J2I J\, the fitting results become 
less satisfactory, presumably because a smaller value of 
rj amplifies effects of both finite system size and higher- 
order terms omitted in the analytic forms (see also the 
discussion below for the estimate of 77). Nevertheless the 
fitting still gives a rather good result at J2/J1 — 1-0 as 
well. This observation gives a strong support to the va- 
lidity of the TLL theory for the SDW2 phase. We empha- 
size that the successful fitting directly demonstrates that 
the characteristic wavenumber of the spin-density wave is 
Q2 = ±7r(l/2 + A7), in accordance with the theory above. 
In the phase diagram in Fig. [Ha), we plot the parameter 
points where the fitting worked well, which cover almost 
the entire region of the SDW2 phase. In the vicinity of 
the 1/3-plateau, however, good fitting results were not 
obtained due to the strong boundary effects. 

In Fig. [3 we present the exponent rj estimated from 
the fitting of the longitudinal-spin fluctuation (sfsf,) — 
(sf)(sf,). Although the estimates have rather large error 
bars coming from high sensitivity to the choice of the data 
range used in the fitting, we can safely conclude that the 
exponent for J2/J1 < 1.5 is always small, i.e., 77 < 0.5. 
This result reflects the fact that the longitudinal-spin 
correlation is the strongest in the SDW2 phase. Fur- 
thermore, the data show the tendency that 77 increases 
with J2/J1. Combining this observation with the per- 
turbative result ([23)1 for J 2 3> J\, we may expect that 
77 increases monotonically with J2/J1 but less than 1 for 
the entire regime of J2/J1 > 1- This means that the 
SDW2 phase, with the dominant longitudinal-spin cor- 
relation, should extend from the intermediate coupling 
regime of J2/J1 ~ 1 to the limit J2/J1 — > 00. 

With decreasing J2/J1, the SDW2 phase appears to 
touch the 1/3-plateau phase. Here we discuss this 
plateau-nonplateau transition within the TLL theory for 
the SDW2 phase. We can consider the effective Hamilto- 
nian with a three-particle umklapp scattering, 

H' + =H+ + \ / rfxsin \ttx{6M - 1) + 3\/2^</) + j , (34) 
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FIG. 7: (Color online) M dependence of the exponent 77 
for the SDW2 phase estimated from the fitting of (sfsfi) — 
i s i){ s i') f° r the antiferromagnetic zigzag ladder with L — 160 
spins. The error bars represent the difference of the estimates 
obtained from the fitting of the data of different ranges. The 
vertical dashed line corresponds to M = 1/6 where the 1/3- 
plateau can appear for small J2/J1. The exponent 77 relates 
to the parameter K+ as 77 = K+ in the TLL theory for the 
SDW 2 phase. 



where %+ is given in Eq. (|27[) , and A is the coupling con- 
stant for three-particle umklapp scattering. The umk- 
lapp term becomes uniform only at M = 1/6. When 
M = 1/6, the umklapp term is relevant for K + < 4/9. 
Then the 0+ field is pinned and acquires a mass gap. 
This results in the 1/3-plateau phase with up- up-down 
spin structure. On the other hand, when approach- 
ing the 1/3-platcau from incommensurate magnetization 
AI -> 1/6, K + takes the universal value K + -> 2/9^ 
This is a commensurate-incommensurate transition. Fig- 
ure [7] indicates that the estimated decay exponent 77 at 
slightly above M — 1/6 seems smaller than 4/9 even 
for J2/J1 = 1.5, suggesting the appearance of the 1/3- 
plateau at this coupling Jij J\. This would mean that 
the upper critical value of the 1/3-plateau phase is larger 
than 1.5, {J2I J\)?2 > 1.5, which is larger than the pre- 
vious estimate (J2/Ji)p2 ~ 1-25 obtained from magne- 
tization curves i2! While our estimated values of 77 may 
contain some large errors, another possible source of this 
discrepancy is that the analysis of magnetization curves 
could miss the plateau with an exponentially small width. 
Further studies with higher accuracy will be needed for 
resolving this issue. 



V. 1/3-PLATEAU PHASE 

The 1/3-plateau phase with a finite spin gap emerges 
at the magnetization M = 1/6 and for the parameter 
regime 0.487 < J2/J1 < 1.25, 21 i 36 i 37 In the 1/3-plateau 



J 2 /J, =b.7,'M=20/120 




FIG. 8: (Color online) (a) Local magnetization (sf) clearly 
shows the up-up-down spin configuration, (b) Semilog plot 
of the absolute values of the averaged correlation functions in 
the antiferromagnetic zigzag ladder with L = 120 spins for 
(J 2 /Ji,M) = (0.7,20/120). 



phase the system has the magnetic LRO of "up-up-down" 
structure , 21 i 36 as shown in Fig. [3Ja). The ground state 
is therefore three-fold degenerate in the thermodynamic 
limit. 

The analysis of magnetization curves has shown that 
the 1/3-plateau phase is surrounded by the TLL1 and 
SDW 2 phases [see Fig. 1 of Ref. M and Fig. E^b) of the 
present paper]. As we discussed in Sees. IIIll and IIV1 we 
can understand this phase diagram as the 1/3-plateau 
phase emerging from instabilities of three-particle umk- 
lapp scattering processes which are inherent in the TLL1 
and SDW2 phases. Here we shall discuss how the up- 
up-down spin configuration emerges through pinning of 
bosonic fields. 

When J2/J1 is small, the plateau emerges from the 
TLL1 phase. As discussed in Sec. lIIH the transition is in- 
duced by the three-particle umklapp scattering process. 
If we fix the magnetization at M = 1/6 and increase 
J2/J1, the umklapp term becomes relevant at K < 2/9. 
Indeed we observed K ~ 2/9 at J 2 /Ji = 0.5 in Fig. [H 
which implies that for J2/J1 > 0.5 the 1/3-platcau phase 
appears. As the umklapp term is relevant, the cf> field is 
pinned at the bottom of the sine potential in Eq. (|20[) . 
v / 47r(0) = tt/2, 7tt/6, and 11tt/6 (A > 0). The bosoniza- 
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tion formula of sf , Eq. ([8]) , then reduces to 



(-l)'asin ( h V4^(^) 



\ 3 



(n = 0,l,2), (35) 



where a > 0. Equation ([35]) gives the up-up-down LRO 
with three-fold degeneracy in the ground state.— ~— 

With larger J%/J\, the plateau phase is next to the 
SDW2 phase. As discussed in Sec. IIV1 this phase transi- 
tion is controlled by the three-particle umklapp term, the 
second term in Eq. (j3"3Jl . When K + < 4/9 and M = 1/6, 
this term becomes relevant, and the <fi+ field is pinned 
to minimize the potential energy. The pinned values are 
V / 2^(0+) = 7r/6, 5tt/6, and 3tt/2 (A < 0). Substituting 
also </>_ = {(f> J) [Eq. ([21))) ] into the bosonized form of sf, 
Eq. H2TJ), yields 



> 1 
s ; = — ha sm 

6 



2ttZ 
IT 



/ 2tt( 



(36) 



which explains the three-fold-degenerate ground state 
with the up- up-down LRO. 

Since both <j) + and </>_ fields are pinned, all low-energy 
excitations in the 1/3-platcau phase are gapped. It thus 
follows that all correlation functions, except the long- 
ranged longitudinal spin correlation, decay exponentially. 
Figure [8] shows the averaged correlation functions for 
J2/J1 = 0.7 and M = 1/6 as a typical example for the 
1/3-plateau phase. The correlation functions decay ex- 
ponentially in accordance with the theory. 



VI. VECTOR CHIRAL PHASE 

The vector chiral phase is characterized by the spon- 
taneous breaking of parity symmetry accompanied by 
nonvanishing expectation value of the vector chirality, 

(«j ) = ((si x si +n ) z ) 7^ 0. The bosonization theory 
for the vector chiral phase was developed in Refs. [32| 
and Hl|, and the appearance of the vector chiral LRO 
in the zigzag spin ladder ([1]) has been numerically con- 
firmed recently 42 ! 43 In this section we present results 
from our detailed numerical study of correlation func- 
tions and compare them with their asymptotic forms de- 
rived from the bosonization theory. 

Let us first briefly summarize the results from the 
bosonization theory. As discussed in Scc. lIVl the effective 
Hamiltonian ([23)) describes the zigzag spin ladder ([1]) in 
the limit J 2 3> Ji- When the g2 term is most relevant, we 
may employ the mean-field decoupling approximation^ 2 , 
in which both d6+/dx and sin(\/27r0_) are assumed to 
acquire nonvanishing expectation values to minimize the 
g<i term. The bosonic fields arc thus pinned as 



where c is a positive constant. Selecting one set of the 
signs from (+,— ) and (— ,+) in Eq. (|3"7) corresponds to 
the spontaneous ^-symmetry breaking in the vector chi- 
ral phase. The antisymmetric sector (</>_, 9—) thus ac- 
quires an energy gap and the low-energy physics of the 
phase is governed by the Gaussian model of the ((/>+, 9+) 
fields, Eq. ([2"T)h in which the 9 + field has been redefined as 
9+ — > 9 + — {d9 + /dx)x to absorb the nonzero expectation 
value of (d0+/dx). The vector chiral phase is described 
by a one-component TLL theory defined by Eqs. (|21[) . 
(ESI, (123), and d37J>. 

Equation (f2"2"|) allows us to write the vector chiral op- 

erators k; 



as 



.(1) 



.(2) 



sin(V27r6>_), 
d9+ 
dx 



(38) 
(39) 



The nonvanishing expectation values in Eq. Q37p result 
in the vector chiral LRO in the ground state. We note 
that the expectation values of the vector chirality satisfy 
the relation 



ID 



2J 2 



„(2) 



0. 



(40) 



so that there is no net spin current^ 2 - Furthermore, one 
can easily obtain the leading asymptotic behaviors of the 
transverse- and longitudinal-spin correlation functions as 
follows: 



= ± 



\r\ 1 / iK + 
A x 



| r |l/4Jf + 



cos(Qr) + • • 
sin(Qr) + 



M 



(41) 
(42) 
(43) 



where Q = (tt+c)/2, and A x is a positive constant. Equa- 
tions (|4Tj) and (|42|) indicate that the spin components 
perpendicular to the applied field have a spiral structure 
with the incommensurate wavenumber Q, which comes 
from the finite expectation value of (d9 + /dx). This he- 
lical quasi-LRO of the transverse components is a char- 
acteristic feature of the vector chiral phase. The sign 
factor ± in Eq. (|42[) comes from the sign ± in Eq. (|37p . 
and it defines the chirality, i.e., the direction of the spi- 
ral pitch. In the longitudinal-spin correlation function, 
the oscillating term with wavenumber Q = 7r(| +M) de- 
cays exponentially as it includes the disordered <j>- field. 
Therefore, if 1/(4K+) < 2, the transverse-spin correla- 
tion function is dominant except the long-ranged vector 
chiral correlations. 

In Fig. [SI we present our DMRG results of the aver- 



aged vector chiral correlation functions 



» K ("') 



for 



{J2I Ji,M) = (1.2,0.35), a representative point in the 
vector chiral phase. Clearly, the vector chiral correlations 
are long-range ordered (the reduction at r > 100 are due 
to boundary effects and should be ignored) . Figure [TU] 
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FIG. 9: (Color online) Averaged vector chiral correla- 
tion functions in the antiferromagnetic zigzag ladder with 
L = 160 spins for (J 2 /Ji,M) = (1.2,0.35). Open cir- 
cles, squares, and triangles respectively represent the DMRG 
data of (fi;Q 1 'ftr 1) }av, -(4 Kr)ay, and (kq 2 ' Kr 2 ')„, where 

(n) / \2 




FIG. 11: (Color online) Averaged transverse-spin correlation 
function (so s r)av in the antiferromagnetic zigzag ladder with 
L = 160 spins for (J 2 /Ji,M) = (1.2,0.35). Open circles 
represent the DMRG data. The fits to Eq. (|41[) are shown by 
solid circles. 



I<K-< 1 V< 1 >> I 
0.03 




FIG. 10: (Color online) Amplitude of the vector chiral cor- 
relations at a distance r = L/2, \{kq K£j 2 )a,v\, in the antifer- 
romagnetic zigzag ladder with L — 160 spins. 




12 J 2/Jl ' A 



FIG. 12: (Color online) J2/J1 dependence of the decay ex- 
ponent 1/(47^+) of the transverse-spin correlation function 
(sosJ:)av in the vector chiral phase for the antiferromagnetic 



zigzag ladder with L = 160 sites. 



shows M and J2/J1 dependences of the amplitude of the 
vector chiral correlations measured at distance r = L/2, 
K^J 2 )a, v \, which indicates the strength of the LRO. 
This figure shows the parameter regions of the vector 
chiral phase; the parameter points where we observe the 
vector chiral LRO are plotted in Fig. [ljb). The vector 
chiral phase appears when J2/J1 is not small, and the 
phase space is split, by the SDW2 phase, into two re- 
gions with either small or large magnetization M. This 
is in contrast with the J\-Ji zigzag ladder with ferromag- 
netic Ji and AF J2 which has the vector chiral phase only 
at small Mi 22 ' 23 ' 53 It is also important to note that each 
one of the vector chiral phases is next to a TLL2 phase 
[see Fig. Hfb)]. The amplitude of the vector chiral order 
parameter exhibits a steep rise at the boundaries to the 
SDW2 and TLL2 phases for small J2 / J\ (see also Fig. 8 
of Ref. IIU for J%/ J\ = 1) while the rise is modest for large 
J2/J1. Incidentally, we have numerically confirmed that 
the vector chiral correlations satisfy the relation (f4"0|). 
These observations on the vector chiral order are consis- 



tent with the previous numerical results i 1 3 

To estimate the TLL parameter K + and the wavenum- 
ber Q of the spiral transverse-spin correlation, we fit the 
DMRG data of (sgs^) av in the systems with L = 120 
and 160 spins to Eq. (jlTj) , with taking Q, K + , and A x 
as fitting parameters. Figure [TT] shows the result for 
(J 2 /J\,M) = (1.2,0.35) and L = 160. We sec that the 
DMRG data are fitted very well to the analytic form, 
except for large distances r > 100 where the boundary 
effect is not negligible. The good agreement between the 
numerical data and the fits supports the validity of the 
TLL theory for the vector chiral phase. 

The decay exponent 1/(4A" + ) of the transverse-spin 
correlation (sgs^) av is shown in Fig. [T^l We have com- 
pared the estimates from L = 160 and 120 spins and 
confirmed that the finite-size effect is negligible in the 
data shown in the figure but not so for some parame- 
ter points (results for which are not shown in Fig. I12[) 
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FIG. 13: (Color online) J2/J1 dependence of the wavenumber 
Q of the transverse-spin correlation function (soSr)av in the 
vector chiral phase for the antiferromagnetic zigzag ladder 
with L — 160 sites. The dashed curve represents the classical 
pitch angle Q = arccos(— J1/U2). 



in the very vicinity of the phase boundaries. It turns 
out that the exponent is rather small, 1/(4K + ) < 1, in 
most parameter region of the vector chiral phase, sug- 
gesting the dominant spiral transverse-spin correlation. 
The exponent becomes larger, as we move closer to the 
1/3-plateau phase. 

Figure Q2] shows the wavenumber Q of the transverse- 
spin correlation function. While Q almost coincides with 
the classical pitch angle arccos(— J1/4J2) near the bound- 
ary to the TLL2 phase, it becomes smaller than the classi- 
cal pitch angle with increasing J2/J1, i.e., moving inside 
the vector chiral phase. We thus find that the incom- 
mensurate wavenumber Q in the vector chiral phase is 
renormalized towards the commensurate value n/2 due 
to quantum fluctuations. 



VII. TLL2 PHASE 

The TLL2 phase is a two-component TLL consisting 
of two flavors of free bosons. In this section, we develop 
its effective low-energy theory based on the bosonization 
of Jordan- Wigner fermions. We then discuss DMRG re- 
sults, which support the effective theory. 

The TLL2 phase is realized in two separated regions of 
high and low magnetic fields in the magnetic phase dia- 
gram. Here we first consider the high- field TLL2 phase, 
for which the origin of the two bosonic modes can be eas- 
ily understood by examining the instability of the fully 
polarized phase. 

Inside the fully polarized phase (h > h s ), the spin- 
wave excitation has a finite energy gap and the dispersion 
relation is given by Eq. (J5|). As the magnetic field is 
lowered, the energy gap decreases and vanishes at the 
saturation field h = h s . For h < h s , the soft magnons 
proliferate and collectively form a TLL. We notice that 
there arc two distinct cases: 



(i) When J2/J1 < 1/4, the bottom of the singlc- 
magnon dispersion is at k = tt (mod 2ir). Magnons with 
k « 7r become soft and condense below the saturation 
field h s , yielding a one-component TLL. Indeed, we have 
found the TLL1 phase in this case (see Sec. IIII[) . 

(ii) When J2/J1 > 1/4, the dispersion has two minima, 
k = it ± Qo with Qo = arccos( J1/4J2). Both magnons 
with k = 7r + Qo and tt — Qq become soft and prolifer- 
ate below the saturation field. The resulting phase is the 
TLL2 phase which consists of equal densities of two fla- 
vors of condensed magnons. We note that, if the densities 
are not equal, the vector chiral phase will be realized,— 
as we will discuss later. 

A similar argument should apply to the TLL2 phase 
appearing at lower magnetic field. The elementary exci- 
tation driving the instability of the dimer ground state 
is a "spinon," a domain wall separating two regions of 
different dimer pattern^-^ For J2/J1 < (J2/Ji)l, the 
dispersion of the two-spinon state has a single minimum 
at k = 7r and only one soft mode is relevant in destabi- 
lizing the dimer state. The TLL1 phase is thus expected 
to show up for M > 0. For J2/J1 > (J2/J1 )l, on the 
other hand, the two-spinon excitation spectrum exhibits 
a double- well structure with minima at incommensurate 
momenta k = ±/co, which leads to the TLL2 (or vector 
chiral) phase for M > 0. The critical coupling at which 
the lowest points deviate from k = n has been estimated 
to be (J 2 /Jr) L - 0.54^ 



A. Two-component TLL theory 

In this subsection we describe the two-component TLL 
theory of the high-field TLL2 phase in detail. As we dis- 
cussed above, this phase can be understood as a two- 
component TLL emerging from condensation of two soft 
magnon modes. This suggests to formulate a low-energy 
effective theory in terms of interacting magnonsi 41 i 76 
Such an approach is valid and useful near the saturation 
field. An alternative approach we adopt here is to for- 
mulate the low-energy theory in terms of Jordan- Wigner 
fermions filling two separate Fermi seas. Advantage of 
the latter approach is that it can be applied in the whole 
TLL2 phase. The connection to the magnon picture will 
also be discussed below. 

We apply the Jordan- Wigner transformation 

s Z l=\-f!fu (44a) 
at = (-1)7/ cxp Uv£/t/J , (44b) 

V n<l / 

*«" = M/// CX P ( 47r E fnfn ) . (44c) 

V n<l / 

to rewrite the Hamiltonian ([I) in the form H = Hq + H' , 
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where 



J2Mj2(f}fl+2 + f} +2 fl 



[2M(J 1 + J 2 )-h]J2flfl: 



and 

//' •/ ^:./;./)::./y. : ./;. : : 

+ ^2 ^ ///; : : fi +2 fi+ 2 : ~fi+ 2 fi : fi+ifi+i '■ 
i 

-flfW-fl+Ji+i-)- (46) 

Here : X : denotes normal ordering of X with respect to 
the filled Fermi sea of fermions with the dispersion 

E(k) = -J 1 cosk+2J 2 M cos{2k)-2M(J 1 +J 2 )+h, (47) 

determined from Eq. (|45[) . Note that the wave number 
k is measured from it as the (—1)' factor is included in 
the Jordan- Wigner transformation. As discussed above, 
in the TLL2 phase the dispersion has two minima and, 
accordingly, there are four Fermi points located at k = 
±k s , ±ki (k s < ki, see Fig. [T4|). The density of fermions 
is 





\ 


1 

/ 


\ 


1 


(45) 


-A 


-kJ 


k\ 


A * 



-(h 

7T 



\-M. 



(48) 



In the limit M — > \, both ki and k s approach Q - In- 
troducing slowly-varying fcrmionic fields for each Fermi 
point, we write the fermion annihilation operator as 



fj = e ik ^ lR (x) +e- ik ' x ip lL (x) 



+ e^ sL (x) + e- ik ^ sR (x) 



(49) 



where the continuous variable x is identified with lat- 
tice index j. We linearize the dispersion around the four 
Fermi points and replace Ho with 

= m J dx (^il^il - ^ir^ir) 

+ 1 v s Jdx 



(50) 



where the velocities vi and v s arc in general different. 
The linearized kinetic term can be written as 



i ~ j 



u—l^s 



dip 



vL 



dx 



( d(fivR 
\ dx 



(51) 



in terms of the chiral bosonic fields ip v ^ and (p uR , which 
obey the commutation relations 



[<p vR (x),ip vR (y)] 
[(p vL {x),<p„ L (y)] 
[(Pv R (x),(p u * L (y)] 



i7rsgn(x - y), 
-i7rsgn(a; - y), 
- -iir6,. 



(52) 



FIG. 14: Four Fermi points located at k — ±k 3 and ±fc;. The 
density of fermions is related to the magnetization, ki — k s = 
7r(i — M). The phase transition to the TLL1 phase occurs 
when k s = 0, i.e., when the inner two Fermi points merge at 
k = 0. 



The fermion densities arc written as 



p uR (x) = :i>t R (x)il> uR (x): = ^^rp 
Pul(x) = :i>l L (x)ip uL (x): = 



(53) 



Finally, the slowly- varying fcrmionic fields are bosonized, 

~>Pu R (x) 



Jb_ e i v „ R (x) 



V2i 



Vvl(z) 



]t__ e -i V „L(x) 



(54) 



where a is a short-distance cutoff on the order of the 
lattice spacing, and r\ v are the Klein factors obeying 
{r]u,riu'} = 25 v y. 

The interaction Hamiltonian H' gives rise to various 
scattering processes of fermionic fields ipuL/R- Among all, 
important in the TLL2 phase are (short-range) density- 
density interactions, 

Hp = it J dx{2g 2l p lL (x)pi R {x) + 2g 2s p sL (x)p sR {x) 

+ 2g 2 j_[piL(x)p sR (x) + piR,{x)p sL (x)] 
+ 9n [piL {x)pil (x) + pir (x)pir (x)] 

+ 9is[p S L{x)p s L{x) + PsR.(x)Psr(x)} 

+ 2g il _[pi L (x)p sL {x) + pir(x)p s r(x)}}, 

(55) 

where g%i s , g 2 ±, 9ai/ s , an d <?4i_ are coupling constants 
that depend on Ji, J 2 , and M. We define the phase 
fields (y = l,s) 

<f) v {x) = -^=[p vL (x) +(p vR (x)], 

(56) 



47T 



Ux) 



1 



VPul{x) - <Pvr(x)]. 



The effective Hamiltonian H 2 = Ho 
quadratic in <f) v and and is diagonalized as 



Hp is then 



(57) 
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by the new fields 8± and <fi± which are linearly related to 
0i s and 4> i a by 



= A> 



Here the 2x2 matrix 



A 



= A~ 



An A w 
A21 A22 



(58) 



(59) 



is a function of the velocities u; s and the coupling con- 
stants <?'s, whose functional form can be found in Ref. [73. 
Without loss of generality, we can assume v+ > V- . 

The Hamiltonian H2 (f57|) is the low-energy effective 
theory of the TLL2 phase. It consists of two free bosonic 
sectors (</> + ,0 + ) and (</>_, 0_). Other interactions which 
arc not included in H p are irrelevant perturbations to 
H2 in the TLL2 phase. An important example of such 
interactions is the backward-scattering interaction 



9i. 



d;r 



^Il{ x )^Ir{ x )^il{x)^ir{x) +H, 

^ / dxcos[V^(9i-9 s )]. 
ir z a z J 



(60) 



The irrelevance of the operator cos[v4tt(0; — 0s)] imposes 
the condition 



1 



(detA) 2 



[(An + A12) 2 + {A 2 i + A22) 2 } >2. (61) 



We note that the vertex operators exp(±iv47T0±) and 
exp(±iv47T0±) have scaling dimension 1. 
The matrix A takes a simple form 



A - J_|V^ 



when the two conditions 







1 1 

1 -1 



vi + gn = v s + g is =: v + 54, 
921 = 92s =: .92, 



(62) 



(63a) 
(63b) 



are satisfied. In this case, the TLL parameters K± and 
the renormalizcd velocities v± are given by 



K± 



V + g 4 ± g4_L ~ g2 T g2J x 

i; + 54 ± 54_l + 52 ± 52J 

2 / _i_ N2-|l/2 



*± = [(« + 34 ± 54±) 2 - (32 ± g2±) 2 ] ' 



(64a) 
(64b) 



This simplified effective theory is applicable when 
J2/J1 >• 1/4 and \M — ^| -C |, i.e., when the magnon 
density is very low and ki — k s -C fc s - In this case one 
can build an effective theory by treating magnons with 
k = 7r±Qo as interacting hard-core bosons i 41 ' 76 We adopt 
a phcnomcnological effective Hamiltonian of interacting 



bosons (0 < v < u)r^ 
H-a = dx 



1 (d%l)\ dip + | dipl dip 
2m 



dx dx dx dx 
+ u{[p+(x)} 2 + Ip-(x)} 2 } + 2vp+(x)p.(x) 

(65) 

where ip_j_ (x) are field operators of two flavors of magnons 
satisfying ft^(x), ip^, (y)] = S^S(x-y), magnon density 

fluctuations p±(x) = ip±(x)?p±(x) — p/2, and m is their 
effective mass. The boson density (per flavor) is assumed 
to be p/2, where p is defined in Eq. (|4"5|). In the low- 
energy, hydrodynamic limit ; 46 ! 79 the magnon fields and 
density fluctuations are written as 



2 



. , 1 dcp±(x) . . 

p±(x) ~ h pcos[TTpx + 2ip±(x)\ 

TT dX 



(66a) 
, (66b) 



where the phase fields obey {tp l _ l (x),d y 'd l _ l i(y)] = 
i^Sfj, :fJ ,'S(x — y). Substituting (|6"6"j) into (p5|) yields 



Hr- 



dx 



[a. \(**± 

\ 4m \ dx 

( df-\ 



+ 



u 
72 



\ dx 



2v dip + dip_ 
tt 2 dx dx 



Once we make the identification of the phase fields, 



dx 

dif- 
dx 



+ vp 2 cos[2(y>+ — </?_)] >. (67) 



f+ = ^{fsL + <Pm), 



IL 



fsR) 



1 



■ <Pm), (68a) 



we can readily see that the Hamiltonian (f6"7| is a special 
case of Hq + H p + H^, with the coupling constants, 



34s = 



9i± 



v 

7T 



vi + gn 

921 = 92 
92± 

Substituting flSHJ} into flM]), we find 



7rp 
4m 



u 

7T 



np 
Am 



u 
n 



v± = x — (u±v), 
" m 



1 ± 1 



/0 



m(u ± u) 



(69a) 
(69b) 
(69c) 

(70) 



Note that v_ and K— vanish when u = v. This 
corresponds to the instability to the vector chiral 
order i 41 ! 76 ' 78 ' 80 We emphasize again that the bosonic ap- 
proach described here is applicable only when | — M <C 1 
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and J 2 /Ji > 1/4, while the general theory (f57|) - (f58|) 
should be valid as a low-energy theory in the whole TLL2 
phase. 

Next we express the spin operators s l using the phase 
fields in the fermionic formulation. We first rewrite the 
string operator used in the Jordan- Wigncr transforma- 
tion, 



\ n<l 



J(k l -k B )x+iy^[fpi(x )+<fi s (x )] 



+ e 



-i(ki-k B )x-i^/Tr[<j>i(x )+c/> a (x )] 



(71) 



where x = x — + , and the second term is added to 
ensure the Hermiticity of the string operator. From Eqs. 
Iggcll , (S3) (EU), and (UU), we obtain 

sj = (-l) x Vs e l ^ 9 ° (x) cos[fc,x + V^Mx)} 

+ (-l) x m e l ^ 6l{x) cos[-k s x + V^Mx)} 
+ (-lfij.e*^' cos[(fc ; - 2k s )x +V^(fa + 20 s )] 
+ (-l) x m e 1 ^ 6 ' cos[(2fci - k s )x +^{2<t>i + fl )] 
+ ..., (72) 

where numerical coefficients are suppressed for simplicity. 
The transverse correlation function becomes 



'0 "r 



— — — cos {kir) + — — — cos (k s r) + 



(73) 

where c/ and c s are constants, and the exponents are 
given by 



x s = 2^12 + ^22) 



1 



(detA) 2 
1 



(dot A) 2 



(74) 



It follows from = ^ — s ; s ; + that 



«f = M-4=^-(0i + 0a) 
V 7T (XX 



+ c\ sin(2fc;a; + v 47T0;) + C2 sin(— 2/c s a; + v 47r</> s ) 
+ c 3 cos[(fc ; - k s )x + Vk{4>i + <j> a )] sm[i/ir(8i - 9 S )] 
+ c 4 cos[(fc ( + k s )x + y/w((f>i - (f>a)] sin[v / 7r(^ - 9 S )} 
+ c 5 cos[2(fc ( - k s )x + \T£k{<\>i + (/>„)] 
+ c 6 cos[2(2fc ; - k s )x + V^r(2(j)i + <f> a )] + ..., (75) 



tained as 



(s z s* r ) = M 2 - 



2it 2 r 2 



[(An + A 12 ) 2 + (A 21 +A 22 )' 



Ci 
\r\ Xl 
C 3 
\r\ Xs 
C 5 



cos(2fc/r) 



C 2 

| r |x 2 



cos(2A; s r) 



cos[(fc; - k s )r] + 
cos[2(fc; — k s )r] ^ 



C 6 



U' ( j 



cos[(fc ( + k s )r] 
cos[2(2fc; — k s )r] 



(76) 



where Cj's are constants, and the exponents are given by 
Xl = 2(A\ l +A 2 l ) 1 
3:3 
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res = 2[(A U + A 12 ) 2 + (A21 + A22) 2 }, 
X(i = 2[(2A n + A 12 f + (2A 21 + A 22 ) 2 ]. 



(77) 



Finally, let us consider local spin polarization (sf ) near 
an open boundary of a semi-infinite spin ladder defined 
on the sites I > 0. Assuming the Dirichlet boundary 
conditions 0/(0) = <j> s (Q) = as in the TLL1 phase [see 
Eq. (|15[)]. we obtain 



(sf)=M- 



ci 



(20 X1 / 2 

C5 



sin(2fc;Z) 



C2 



(20 X2 / 2 



sin(2fc s Z) 



+ 



(2iys/ 2 
(2iy*/ 2 



cos[2(fc/ — k s )l] 
cos[2(2fc; - k s )l] 



(78) 



Observe that the exponents in Eq. ([75)1 are a half of the 
corresponding ones in Eq. (|76[) and that the vertex oper- 
ators of the 9 1 yS fields do not contribute to Eq. (|75|) . 

An important characteristic feature of the spin corre- 
lations ([73")) and (|T6"|) in the TLL2 phase is the presence 
of two incommensurate (Fermi) wave numbers ki and k s 
(and their linear combinations). 

Before closing this subsection, we note that Frahm and 
Rodenbeck studied an exactly solvable zigzag spin ladder 
model with additional three-spin interactions.—^ Their 
model has a phase corresponding to our TLL2 phase. 
They have calculated, using the Bethe ansatz solution 
and conformal field theory, exponents of several terms in 
the longitudinal spin correlation (|76l) . 



Instabilities 



where Cj's are nonuniversal constants. The long-distance 
behavior of the longitudinal spin correlation is then ob- 



In the magnetic phase diagram (Fig. [I) each TLL2 
phase is next to a vector chiral phase and the 
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TLL1 phase. Since these neighboring phases are one- 
component TLLs, one of the two massless modes in the 
low-energy Hamiltonian (f57|) has to become massive or 
disappear from low-energy spectra at the transitions from 
the TLL2 phase. Here we discuss instabilities of gapless 
modes in the TLL2 phases which cause the phase transi- 
tions to the vector chiral and TLL1 phases. 

As pointed out by Kolczhuk and Vckuar^ in the inter- 
acting magnon picture valid in the vicinity of the satu- 
ration field [Eqs. (p5)) -(fT0" |) ], the instability to the vector 
chiral phase corresponds to the "demixing" or "phase 
separation" instability ; 78 ' 80 which occurs when both V- 
and K- vanish. Alternatively, if we regard the two fla- 
vors as up and down pseudospins, the TLL2 and vector 
chiral phases correspond to paramagnetic and ferromag- 
netic phases, respectively. The transition between the 
TLL2 and vector chiral phases is then regarded as a fer- 
romagnetic transition^ Away from the saturation field, 
the interacting magnon picture is no longer applicable, 
and we should use the low-energy effective Hamiltonian 
(f57|) with the A matrix (f59|) . The instability to the vector 
chiral phase is then signaled by v_ — and det A = 0. 

The transition between the TLL2 and TLL1 phases 
is characterized by a cusp singularity in the magnetiza- 
tion curved Since M and h correspond to the particle 
density and the chemical potential of the Jordan- Wigner 
fermions, the origin of the cusp singularity can be at- 
tributed to the van Hove singularity of the fermion den- 
sity of states, which exists at the saddle point k = of the 
dispersion (|4"T|) . Thus, the TLL2-TLL1 transition is con- 
sidered to occur when the chemical potential matches the 
saddle-point energy, and the two Fermi seas merge into 
a single Fermi sea<^ Indeed, the Bethe-ansatz study of a 
solvable model finds that the transition of commensurate- 
incommensurate type occurs when k s = 0i 81 ' 82 In our 
low-energy effective theory, the transition is driven by 
the operator (the C2 term in sf), 

hj dxsm(-2k s x + V^ncps), (79) 

which turns into a mass term (scaling dimension 1) for 
fermions at the TLL2-TLL1 transition. Comparison 
of our effective theory with the Bethe-ansatz study in 
Ref. HH shows that the A matrix takes the form 

at the transition (h \, h c {), in agreement with our pic- 
ture of the TLL2-TLL1 transition as a commensurate- 
incommensurate transition caused by the operator ([79| . 
Here £ is the dressed charge defined in Ref. |82| . 

C. Numerical results 

In Fig. [2^e) we have shown the correlation functions 
at J2/J1 = 0.6 and M = 0.4, as a typical example of 



the TLL2 phase. We see that both the longitudinal- and 
transverse-spin correlation functions decay algebraically. 
The vector chiral LRO is clearly absent. 

As we have discussed in Sec. IVII A[ the defining fea- 
ture of the TLL2 phase is that its low-energy physics is 
governed by the two independent sets of free bosons. The 
low-energy theory is a conformal field theory with central 
charge c = 1 + 1. The central charge can be numerically 
measured through the entanglement entropy, 

S(l) = -Tr a [p(l)lnp(l)], (81) 

where the reduced density matrix for the subsystem fi = 
{ s j \ 1 < j < is defined by 

p(l) =Tr n \0)(0\. (82) 

Here |0) is the ground state wave function, and the spins 
sj+i, . . . , sl in the environment f2 are traced out. The 
entanglement entropy of a ID critical system with open 
boundaries is known to have a logarithmic dependence 
on i^2£ 

S(l) = -In/ + const., (83) 
6 

in the thermodynamic limit, L — > 00 and / 1. For 
finite-size systems of L spins, Inl in Eq. (|83[) should be 
replaced by^ 



x = In 


~ L • 1 




— SIB 






7T ' 





(84) 



Hence we can measure the central charge c as a coefficient 
of x. This method was recently used to detect the central 
charge of the critical spin Bosc metal phase in a related 
model of the J\-Ji zigzag ladder with a ring exchange 
interaction^ 7 - Figure Eta) shows the entanglement en- 
tropy S(l) in the TLL2 phase (J2/J1 = 0.6, M = 0.4) 
as a function of x. We clearly see that S(l) ~ x/3, in- 
dicating that c = 2. For comparison, we have computed 
the entanglement entropy in the TLL1 and vector chi- 
ral phases. The numerical results shown in Fig. ITST b) 
demonstrate that S(l) ~ x/6 for large x, i.e., c = 1. 

Having confirmed that the TLL2 phase has c = 2, i.e., 
that the low-energy physics is governed by two free bo- 
son theories, we now discuss spin correlation functions. It 
turned out, however, that the presence of the two Fermi 
wavenumbers ki and k s makes it difficult to analyze cor- 
relation functions in the TLL2 phase. For this reason we 
focus attention to the simplest, one-point function (sf). 
The Friedel oscillations near open boundaries give us in- 
formation on the Fermi wavenumbers. 

We show in Figs. [H] and [T7] the squared modulus of 
the Fourier transform of the local spin polarization 

*"(*) = -7fE e< *'« fl ">- M )- ( 85 ) 

At J2/J1 = 0.6 (Fig. [16]) the TLL2 phases appear when 
< M < 0.075 and 0.25 < M < 1/2, and the TLL1 
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FIG. 15: (Color online) (a) Entanglement entropy of the 
TLL2 phase (J2/J1 = 0.6, M = 0.4). The horizontal axis 
is x = ln[(L/7r) sin(7rZ/L)]. Solid and dotted lines represent 
the slope 1/3 (c = 2) and 1/6 (c = 1), respectively, (b) 
Entanglement entropy of the TLL1 phase (J2/J1 = 0.1, M = 
0.375) and the vector chiral phase (J2/J1 = 1.2, M = 0.35). 
Dotted lines indicate the slope 1/6 (c = 1). 
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FIG. 16: (Color online) (a) Squared modulus of the Fourier 
transform of the local spin polarization, |s z (fc)| 2 , for the 
antiferromagnetic zigzag ladder with L — 160 spins and 
J2/J1 = 0.6. (b) M dependence of peak positions of \s z (k)\ 2 . 
Solid line represents the highest peak while the dotted lines 
correspond to the subdominant peaks in TLL2 phase. Gray 
horizontal lines show phase boundaries. 



phase is located at 0.1 < M < 0.2. In the TLL1 phase 
we see a very sharp peak in |s z (fc)| 2 at k — 7r(l — 2M), 
in agreement with Eq. f|15[) . Although greatly reduced in 
magnitude, the peak persists in the TLL2 phases. This 
faint peak comes from the fourth term, with wavenum- 
ber 2(ki — k s ), in Eq. (|78]l. We attribute the strongest 
peak of |s z (/c)| 2 in the TLL2 phase to the second term 
in Eq. (|78[) with wavenumber 2/q. The two peaks meet 
when the TLL2 phase is turned into the TLL1 phase, 
i.e., when k s vanishes, in accordance with the discus- 
sion in Sec. I VII Bl Moreover, at the saturation limit 
M — > 1/2, the wavenumber k max of the strongest peak 
approaches 2Q = 2 arccos( Ji/4J 2 ), where Qo is the mo- 
mentum of the soft magnon in the fully polarized state, 
while /c max — > 2fco as M — > 0, where ko is the mo- 
mentum of the soft singlc-spinon excitation in the dimcr 
phase estimated numerically^ In the higher-field TLL2 
phase we see a third, faint peak, whose wavenumber k^ 
equals 2Q at M — > 1/2 and increases with decreasing 
M . We have found numerically that k% — k max equals 
2(ki — k s ) — 7r(l — 2M) modulo 2-7T, from which we con- 
clude /C3 = Aki — 2k s . Interestingly, |s z (A:)| 2 does not 
have a peak corresponding to k = 2k s . Comparing the 



peak heights, we can deduce the following inequalities for 
exponents, 

Xl < X 5l X 6 < X 2 ,X 3 ,X4. (86) 

From the relation xjx s = xi/x 2 , we can also obtain 

xi < x s . (87) 

These observations suggest that the dominant component 
in the transverse-spin correlation function comes from the 
first term in Eq. (|73j) with a wavenumber ir±ki while the 
dominant longitudinal-spin correlation comes from the 
third term in Eq. (|T6[) with a wavenumber 2fc;. 

Figure [T7l shows |s z (fc)| 2 at J2/J1 = 0.9. In this case 
we have the TLL2 phase for 0.35 < M < 1/2, the SDW 2 
phase for 0.075 < M < 0.2, and the vector chiral phase 
for < M < 0.05 and 0.2 < M < 0.3. Characteristics 
of incommensurate wavenumbers giving rise to the peaks 
in |s z (fc)| 2 in the TLL2 phase are the same as in Fig. [16] 
In the SDW2 phase the strong peak is found to be at 
k = tt(1/2 + M), in agreement with Eq. 
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FIG. 17: (Color online) (a) Squared modulus of the Fourier 
transform of the local spin polarization, |s z (A:)| 2 , for the 
antiferromagnetic zigzag ladder with L — 160 spins and 
J2/J1 = 0.9. (b) M dependence of peak positions of |s z (fc)| 2 . 
Solid line represents the highest peak while the dotted lines 
correspond to the subdominant peaks in TLL2 phase. Gray 
horizontal lines show phase boundaries. 



model is critical. While the interladdcr couplings can 
have a complicated geometry, it is quite natural to ex- 
pect, to the first approximation, that the ladders are 
coupled in a non-frustrated way. In such a case, the dom- 
inant algebraic correlation in the purely ID model leads 
to the magnetic LRO in the real quasi-lD compounds. 
Based on our results on the correlation functions, we 
can thus predict that several different magnetic-ordered 
phases appear in the quasi-lD zigzag ladder compounds; 
In the parameter regime of the TLL1 phase, we expect a 
canted antiferromagnetic ordered phase for small J2/J1 
and an incommensurate longitudinal spin-density wave 
ordered phase with a wavenumber Q = 7t(l ± 2M) for 
slightly larger Ja/Ji- The region of the SDW2 phase 
will be replaced by an incommensurate longitudinal spin- 
density wave ordered phase with Q2 = ±7r(l/2 + M). 
The vector chiral phase turns into the spiral ordered 
phase, in which spins perpendicular to the applied field 
have incommensurate long-range order. This is similar to 
the classical helical magnetic structure albeit with renor- 
malizcd pitch and canting angles. For the parameter 
regime of the TLL2 phase, the system should exhibit 
the coplanar "fan" phase characterized by the coexis- 
tence of incommensurate longitudinal- and transverse- 
spin LROs. This is consistent with the argument by 
Ueda and Totsuka^ they showed, using a dilute Bose 
gas description, that the coplanar fan phase appears near 
saturation in the quasi-lD system in a wide parameter 
region around J2/J1 — 1/3. 



VIII. CONCLUDING REMARKS 

By the thorough comparison between numerically ob- 
tained correlation functions and asymptotic behaviors 
derived from low-energy effective theories, we have iden- 
tified the nature of critical TLL phases that appear in 
the spin-1/2 J1-J2 AF Hcisenbcrg zigzag ladder under 
magnetic field. These critical phases consist of three 
one-component TLL phases (the TLL1, SDW2, and vec- 
tor chiral phases) and a two-component TLL phase, the 
TLL2 phase. From the fitting, we numerically estimated 
the TLL parameter in one-component TLL phases as a 
function of J2/J1 and the magnetization A/. The results 
allow us to determine the decay exponents of the alge- 
braic spin correlation functions and reveal the dominant 
correlation function in each phase. In addition, we de- 
veloped an effective theory for the two-component TLL, 
which reasonably reproduces numerically obtained cor- 
relation functions in the TLL2 phase, which appears in 
two parameter regions in between the TLL1 and vector 
chiral phases. 

One of important implications of our results concerns 
field-induced phase transitions in quasi-lD compounds, 
in which weak interladder couplings usually induce a 
magnetic LRO when the ground state of the pure ID 



Another related quasi-lD system is a spatially 
anisotropic triangular antiferromagnet, with interchain 
exchange J' much weaker than the intrachain exchange 
J. This model was studied recently*^ and the obtained 
phase diagram shows a resemblance to that of the zigzag 
ladder. In ID limit of J' -C J, Starykh and Balcnts 9 
found a collincar spin-density wave with wave vector 
k x = 7r(l ± 2M) in intermediate magnetic field regime 
and a cone phase with spiral transverse order in high 
magnetic field regime. Kohnoii also found instability to 
the ordering of incommensurate longitudinal spin-density 
wave with momentum k x = 7t(l ± 2M) applying weak- 
coupling analysis to ID exact solution. If we take a zigzag 
ladder out of this anisotropic triangular system, the na- 
ture of the incommensurate spin-density wave and cone 
phases, respectively, is essentially the same as that of the 
SDW2 and vector chiral phases we showed in the regime 
of Ji <C Ji- (Note that the definition of the unit length 
along chains on the anisotropic triangular lattice is twice 
larger than that we used in the zigzag ladder.) Tran- 
sitions from the cone phase to coplanar fan phase with 
increasing J'/ J were also discussed in Ref.|H, which pre- 
sumably relate to the transitions from the vector chiral 
phase to the TLL2 phase with increasing J1/J2 in the 
zigzag ladder. 
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